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Abstract 

We extend to the strange quarks and antiquarks, the statistical approach of 
parton distributions and we calculate the strange quark asymmetry s — s. We 
find that the asymmetry is small, positive in the low x region and negative 
in the high x region. In this framework, the polarized strange quarks and 
antiquarks distributions, which are obtained simultaneously, are found to be 
both negative for all x values. 
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1 Introduction 



Although strange quarks and antiquarks s and s play a fundamental role in 
the nucleon structure, they are much less known than the parton distribution 
functions (PDF) for the light quarks u and d. The measurements of the 
structure functions in deep inelastic scattering (DIS) of charged leptons on 
hadrons provide the best informations on u, d, whereas neutrino DIS and 
lepton-pair production in hadron collisions put some constraints on the sea 
quark distributions u and d. Concerning the strange quarks, due to the fact 
that the structure functions are largely dominated by u and d, the extraction 
of the small components s and s is rather difficult. Therefore most of the 
phenomelogical models for the PDF studies use the simplifying assumption 
s{x) = s{x) = k{u + d)/2 (with k ~ 0.5). However, nothing prevents s(x) ^ 
s{x) and we will see how to achieve this inequality in the statistical parton 
model [H [21 111]. 

An experiment on neutrino (antineutrino) -nucleon charged-current DIS 
by the CCFR collaboration [5] at the Fermilab Tevatron has measured the 
production of dimuon final states coming from a charm quark fragmentation. 
This process involves the interaction of a neutrino (antineutrino) with an s 
(s) or d {d) quarks, via a exchange, which can be used to isolate their 
distributions. Since the contribution of the down to charm production is 
Cabibbo suppressed, scattering off a strange quark is responsible for most 
of the total dimuon rate. Unfortunately, only an average value of s + s was 
extracted from this experiment, but the size of strange quark distribution was 
known for the first time. Later, the NuTeV collaboration [6] has reached a 
greater accuracy by a high-statistics dimuon measurement, allowing to study 
independent information on s and s and the difference s — s. 

On the theoretical side one of first attempt to separate the s and s distri- 
butions was investigated in a light-cone model |7] and more recently, other 
models based on nonperturbative mechanisms were proposed [H [9] . A global 
QCD fit to the CCFR and NuTeV dimuon data has shown a clear evidence 
that s s [lU]. In another approach based on perturbative evolution in QCD 
at three loops [TT], one is able to generate a strange-antistrange asymmetry 
although at the input scale s = s. 

In this Letter, we will show how to construct the strange and antistrange 
quark distributions in the statistical parton model. Since according to our 
method, the basic distributions are the helicity dependent ones, and s^, 
we will obtain simultaneously the unpolarized, s = + s~ , s = + s~ , 
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and the polarized PDF, As = — s~ , As = — s~ . We will also explain 
how to determine the few parameters involved. Our results will be compared 
with other theoretical predictions. 

2 Strange quark and antiquark distributions 

In the statistical approach the nucleon is viewed as a gas of massless partons 
(quarks, antiquarks, gluons) in equilibrium at a given temperature in a finite 
size volume. Like in our earlier work on the subject [1], we propose to 
use a simple description of the parton distributions p{x) proportional to 
[exp[(a; — Xop)/a;] ± l]""*^, the plus sign for quarks and antiquarks, corresponds 
to a Fermi-Dirac distribution and the minus sign for gluons, corresponds 
to a Bose-Einstein distribution. Here Xop is a constant which plays the 
role of the thermodynamical potential of the parton p and x is the universal 
temperature, which is the same for all partons. Since quarks carry a spin-1/2, 
it is natural to consider that the basic distributions are qf{x), corresponding 
to a quark of flavor i and helicity parallel or antiparallel to the nucleon 
helicity. From the chiral structure of QCD, we have two important properties 
which allow to relate quark and antiquark distributions and to restrict the 
gluon distribution: 

- The potential of a quark of helicity h is opposite to the potential of the 
corresponding antiquark of helicity -/i, therefore X^^ = —X^^. 

- The potential of the gluon G is zero Xqq = 0. 

The sum rules, coming from the quantum numbers of the proton, u — u = 2 
and d — d = 1, give rise to higher values for the potentials of the u's than 
for the d's. In fact we have found Xq^ > Xq^ ~ Xq^ > X^^, which is 
also consistent with the known facts that Au{x) > and Ad{x) < 0. This 
ordering leads immediately to some important consequences for the light 
antiquarks, namely 

i) d{x) > u{x), the flavor symmetry breaking, which also follows from 
the Pauli exclusion principle, whose effects are incorporated in the statistical 
model. 

ii) Au{x) > and Ad{x) < 0. 

We now turn to the procedure to construct the strange quark distribu- 
tions. In the original version of the statistical parton model [1] we have 
assumed that the unpolarized strange quark and antiquark distributions are 
equal and they can be described by a linear combination of light antiquark 
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distributions at the input scale Ql, namely 

xs{x, Ql) = xs{x, Ql) = ^ [xu{x, Ql) + xd{x, Ql)] , (1) 

where the coefficient 1/4 is an average value of some current estimates. For 
the corresponding polarized distributions a similar assumption was made, 
more precisely 

xAs{x, Ql) = xAs{x, Ql) = ^[xAd{x, Ql) - xAu{x, Ql)] , (2) 

which leads to a large negative distribution, since Ad < and Au > 
(See Fig. 18 of Ref. p]). In order to introduce a difference between s and 
s, here we follow the procedure used earlier to built the light quarks PDF, 
with the recent improvement obtained from the extension to the transverse 
momentum of the PDF [4]. So the strange quark distributions s^{x,Ql) of 
helicity h = ±, at the input energy scale Ql = 4GeV^, have the following 
expressions 

, ^ AX+y^ ln(l + exp [kX^Jx]) 

^ exp[(x-Xo^J/x] + l ln(l+exp[A;Xo+/x]) exp(x/x) + 1 ' 

n 

and similarly □ for the antiquarks s^{x, Ql) 



b 



exp[(x + Xq-/) /x] + 1 In (1 + exp [-kX^Jx]) exp(a;/a;) + 1 ' 

(4) 

The value of the input energy scale is arbitrary and should not affect the 
results which satisfy the Q^ QCD evolution. Our choice was dictated in 
Ref. [1] by the existence of many accurate data at Q^ = 4GeV^. The ffist term 
in the right hand side corresponds to the non diffractive part, while the second 
is associated with a diffractive component common to all distributions. The 
ratio of the logarithms originates simply from our extension of the statistical 
distributions to the transverse degree of freedom and justifies the presence 
of a multiplicative factor in the Fermi-Dirac functions, ffist introduced in 
Ref. |1J, as was explained in Ref.|l]. The above expressions involve some 



^As mentioned above, quarks and antiquarks are not independent due to the relation 
between the potentials Xq^ = —X^^. 
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parameters already determined in our previous works [H H] , which we recall 
now 

A = 1.74938, A = 1.90801, Xo+ = 0.46128, X^^ = 0.22775, 

X = 0.09907, b = -0.25347, k = 1.42 . (5) 

Therefore it remains only four free parameters to determine, namely the 
two potentials X^^, bs and the normalization of the diffractive part Ag. 

In order to obtain these free parameters, we will use some constraints: 
first, the nucleon does not have any strangeness quantum number, as a con- 
sequence the asymmetry [a~] has to vanish for all x values 

[a~] = / [s{x) — s{x)]dx = , (6) 
Jo 

second, from the second Bjorken sum rule, the first moments of the polarized 
quark distributions must satisfy the relation 

Aqs = Au + Au + Ad + Ad- 2{As + As) = 3F - D , (7) 

where F and D are the hyperon beta decay constants, so that 3F — D = 
0.579 ± 0.008. From the values of the first moments of the light quarks 
calculated in Ref . [T] , where quarks and antiquarks are related through their 
potentials, we can deduce another constraint, namely 

[a+] = [ [As{x) + As{x)]dx = -0.04675 . (8) 
Jo 

From the above discussion on the light quarks, it is now clear that the sum 
rule Eq. will lead to strange potentials X^^ smaller than Xq^ and X^^. 
This obvious expectation has been observed in several earlier works on the 
same subject, see for example Ref. [12]. Similarly in order to satisfy Eq.(8), 
we anticipate that we will find Xq^ > X^^. To determine the free parameters, 
in addition to the above constraints, we will use some experimental results 
obtained by the CCFR and NuTeV collaborations [SI E] on the production 
of dimuons from neutrino and antineutrino scattering on iron. We have 
performed a next-to-leading order (NLO) QCD analysis of the data, keeping 
the light quark distributions as in Ref. [Tj. We have obtained the following 
values for the parameters, Xo+ = 0.08101, Xg" = 0.20290, 6, = 2.05305 and 
As = 0.05762, all with an error of the order of few percents. We observe that 
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the chemical potentials for the strange quarks are smaller than the potentials 
for light quarks u and d and that Xq^ > X^^, like in the case of the d quark. 
Due to the large value of bg, the contribution of the non diffractive component 
is strongly suppressed, in the small x region. 

Our results are displayed in Figs. 1-5 and the fit is rather satisfactory since, 
as an indication, we have a x^/dof of the order of 1.5, compared to 2.75 if 
one uses instead, the simplifying assumption Eq. (1). We have also checked 
that in this ealier version, it is not possible to reproduce the rapid rise of 
s + s at low X and = 4GeV^ of the data, as shown in Fig. 1. 

At the input scale Qq = 4GeV^, [a~] = is satisfied to a great accuracy 
and we have checked that this constraint is not affected by the evolution. 
We also find [a"*"] = —0.0221, which is compatible with a recent HERMES 
determination [T3], for this first moment in the measured region, x > 0.02 
and < >= 2.5GeV^ namely 0.006 ± 0.029(stat.) ± 0.007(sys.). 

For the first moment of the asymmetry 



we have [S-] = -0.00194, to be compared with the value -0.0011 ± 0.0014 
found by the NuTeV collaboration [H] , and with the allowed range extracted 
from a global QCD fit by CTEQ [lOj -0.001 < [S'] < 0.004. The calcu- 
lations in the light-cone meson-baryon model, lead to two positive results, 
namely 0.0042 < [S~] < 0.0106 for the choice of a Gaussian wave function 
or 0.0035 < [S~] < 0.0087 for a power-law wave function [9J. 

We now turn to discuss the main features of the distributions obtained 
from this fit and compare them with other theoretical models. We show 
the unpolarized and polarized strange quark distributions at the input scale 
Qq = 4GeV^ in Fig. 6. We observe that the distributions s{x) and s{x) are 
almost identical for x < 0.05, because the diffractive component dominates 
largely, whereas s{x) is a little larger than s{x) for 0.05 < x < 0.25 and 
s{x) < s{x) for 0.25 < x < 1. These features remain unchanged for higher 

values, as shown in Fig. 7 for the difference s — s plotted as a function 
of X, for = 4, 20, lOOGeV^. This pattern is similar to that one gets in the 
meson cloud model [8jJ and also in another approach based on perturbative 
evolution in QCD at three loops [11], although, in this later case the sign 
change occurs at a much smaller value of x. On the contrary, in Ref.[9] they 
found that s{x) < s{x) in the small x region and s{x) > s{x) in the large x 
region. 
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Finally, both As{x) and As{x) are negative for all x, as shown in Fig. 6, 
in reasonable agreement with the results of Ref. [12]. This contradicts the 
expectation of the meson cloud model [8], so it is clear that we need for a 
better measurement of the strange quark contribution to the nucleon spin, 
has was also stressed in Ref. [T5] . 

3 Conclusion 

We have investigated the possibility to introduce an asymmetry for the 
strange quark distributions in the framework of a statistical parton model. 
In the absence of direct precise experimental data, we have imposed different 
unpolarized and polarized constraints on the distributions and an extensive 
use of the recent results from CCFR and NuTeV. The main results are that 
s{x) — s{x) is indeed small, as expected, positive in the low x region and 
negative for x > 0.25. Our approach has the unique feature to provide si- 
multaneously the polarized distributions for strange quarks and antiquarks 
which are found to be both negative for all x. New results on the strange 
quarks distributions are welcome, because they will produce further tests on 
the present determination and hopefully some improvement on them. 
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Figure 1: The average strange quark antiquark distributions determined at 
NLO as a function of x for = 4^ 20, lOOGeVl Data from CCFR Collabo- 
ration [5]. 
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Figure 2: Comparison of the CCFR v data [6] to the result of the fit for 
da/dxdy, in units of charged- current a, for various kinematic ranges in en- 
ergy, X and y. 
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Figure 3: Comparison of the CCFR v data [6] to the result of the fit for 
da jdxdy in units of charged-current a, for various kinematic ranges in energy, 
X and y. 



11 




Figure 4: Comparison of the NuTeV v data [6] to the result of the fit for 
da / dxdy , in units of charged- current cr, for various kinematic ranges in en- 
ergy, X and y. 
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Figure 5: Comparison of the NuTeV v data [6] to the result of the fit for 
da / dxdy , in units of charged- current cr, for various kinematic ranges in en- 
ergy, X and y. 
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Figure 6: The unpolarized and polarized strange quark and antiquark distri- 
butions determined at NLO as a function of x for = 4GeV^. 
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